Josephson vr-state in superconductor-Luttinger liquid hybrid systems 
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Josephson current through a Luttinger liquid (LL) under a magnetic field is theoretically studied. 
We derive an analytical expression of Josephson current for clean interfaces, by using quasiclassical 
Green's function and functional bosonization procedure. We show that critical currents can be 
renormalized by electron-electron interactions at perfect transparency when LL is adiabatically 
connected with superconductors. We also find that a generation of 7r-state, due to spin-dependent 
energy shift in Andreev bound states (ABS), is prohibited even at zero temperature when the 
strength of repulsive interactions reaches some critical value. The suppression of vr-state is caused 
by the low energy fluctuations propagating in LL, and making the Zeeman splitting in ABS blurred. 

PACS numbers: 71.10.Pm, 71.70.Ej, 74.50.-|-r 
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Transport phenomena in one-dimensional (ID) struc- 
tures are strongly affected by electron-electron interac- 
tions. Fermi liquid description breaks down due to elec- 
tron correlation, and systems are believed to behave as 
Luttinger liquids (LLs), where low- lying excitations are 
collective modes rather than single-particle excitations. 
This gives rise to characteristic phenomena such as spin- 
charge separation and charge fractionalization. These 
behaviors have been verified experimentally, e.g., in car- 
bon nanotubes (CNTs) P, Q and quantum wires on 
GaAs/AlGaAs heterostructure 0, 0. Inspired by ad- 
vances in microfabrication techniques, hybrid structures 
connecting LL with other conductors have attracted at- 
tentions from theoretical and experimental sides as can- 
didates of new mesoscopic electronic devices. 

Recently, superconducting proximity effect on CNTs 
in contact with superconductors was reported in a cou- 
ple of experiments |^ 0| . There are some previous works 
about Josephson currents in superconductor-Luttinger 
liquid (S-LL) hybrid systems ill 1^, ll ES IH Dl] • In such 
systems, electron-electron interactions modify the nor- 
mal and the Andreev reflection (AR) process, and lead 
to unique power law behaviors of low energy properties 
such as proximity effect, conductance and local density 
of states (LDOS). On the other hand, a generation of 
metastable 7r-state has been studied theoretically ^3] and 
experimentally |l4l Il5| in superconductor-ferromagnct 
(S-F) junctions. The 7r-state is a result of spin-dependent 
energy shift in Andreev bound states (ABSs) due to ex- 
change energy in F. Here, we study Josephson current in 
the presence of the 7r-state in S-LL junctions. By apply- 
ing a magnetic field to the LL region, we can consider 
that Zeeman energy in LL plays the same role as the ex- 
change energy in F. In this paper, we derive the analytical 
expression of Josephson current and find that critical cur- 
rents can be renormalized by the interactions at perfect 
transparency when LL is adiabatically connected with 
reservoirs. We also show that the Coulomb interactions 
affect the generation of the metastable 7r-state through 
the LDOS in ID structures. 

The system under consideration is a long S-LL-S junc- 
tion where LL is adiabatically connected with s-wave su- 
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FIG. 1: A possible implementation of system setup: s-wave 
superconductors are deposited on quantum wire which is adi- 
abatically connected with two dimensional Fermi liquids. We 
believe that superconducting proximity effect penetrates in 
Fermi liquid regions. The electrons in superconductors can 
tunnel into Luttinger liquid via Andreev reflection process. 
For sake of simplicity, barrier potentials at the junctions are 
assumed to be negligible. 



perconductors (see Fig.l). Approximately influence of 
the interactions in superconductors is neglected, and su- 
perconducting order parameters IS.e^^^^'^ is assumed to 
change abruptly. For LL region, we assume that Fermi 
velocity is spin-independent vp] = vpi = vp- This as- 
sumption may be justified as far as the Zeeman energy is 
far smaller than the Fermi energy. In this case, we can 
consider that the spin-charge separation in LL does not 
break down [T^ . Thus, Hamiltonian with usual g-ology 
is expressed as 



i?LL = / dx 



d 

+ X! ^/'i,s(a;,t)pa,s'(a;,i) 



(1) 



where pa^s and 
ergy. Here, a = 



h are particle density and Zeeman en- 
± and s = ± denote direction of move- 
ment and spin, respectively, p = {||,_L} specifies in- 
teraction between electrons with parallel or anti-parallel 
spins. We consider the case where back-scattering and 
umklapp-scattering processes are irrelevant for simplic- 
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ity. Throughout this paper, h and fc^ are set to unity. 

We apply a so-cahed Hubbard-Stratonovich transfor- 
mation to the interaction terms ^3 , and obtain the ac- 
tion using auxihary field <j)a^s{x,t) as follows 



S[4>]^ j dtdxy Loiipl^p) +Li{(t>) 

+ ^<l>a.,s{x,t)pa.s{x,t) . (2) 

a,s 

Here, Lq and Li are the Lagrangian densities of free 
fermions and collective fluctuations; they can be ex- 
pressed in matrix form as 

9 ^ d - 

Lo = [ffj^i'i-Q^'^ + ^'"F-Q^Tz^z + hT^)^lj], (3) 

where g represents interactions, and 



(4) 
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with (Ti being usual Pauli matrices. Throughout the pa- 
per, quantities with gcareth denotes (4 x 4) matrices, and 
those with boldface (2 x 2) matrices. 1st and 3rd row cor- 
respond to right and left moving electrons with spin up, 
and 2nd and 4th row to left and right moving holes with 
spin down. From Eq.(2), one can regard LL as the free 
fermions propagating in the bosonic "environment" [Tsj . 
Here the free fermion part includes the role of topological 
term in usual bosonization method. Therefore, we treat 
(j){x, t) as "local scalar potential" for a time, and average 
it in terms of fluctuation fields. One can obtain quasi- 
classical Green's function in LL by solving Eilenberger 
equation ^ 



d 

ivF-g^g{x,t,t'\(j)) + 

d 

{i-g^ + h)T2T,z + (j){x,t)Y,z, g{x,t,t'\ 



= 0,(6) 



where [■■■]- denotes commutator. 

In a similar fashion, one can construct the stationary 
Green's functions in S regions Here, it is assumed 
that the magnetic field is not applied to S region. In 
Fourier representation with respect of time difference, 
they are given by 



g^^^\x,e) = U{go{e)+g'^^''Hx,e))Ul 



(7) 



where U = exp[ixTz/4]- 3o(e) = (e'Tz -I- iAfj^)/f2(e) is the 
one at far from the interfaces, and 
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1 a(e) 
-a(e)Ti -1 



±1 I , (8) 



where L (R) depicts S of left (right) hand side. Here, 
n{e) = v/e^ - A2 and a(e) = A/{e + n{e)). One can 
see that a{e) is equal to the AR amplitude for inter- 
faces with non-interacting normal conductor. We deter- 
mine coefficients C^^^\e) by applying boundary condi- 
tions at the interfaces x = ±i/2. Since we focus on 
the junction with clean interfaces, the boundary condi- 
tions are that they are continuous along a trajectory, i.e., 
5(±f - 0, i, t'|0) = .g(±f + 0, t, t'\ct)) p. 

If L is larger than coherence length f = vf/ttA, we 
find that ABS has the following form 
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where n is an integer, and 7 = 2Lh/vF (mod 27r) is phase 
shift due to the magnetic field. Here, the influence of the 
auxiliary field appears as a small energy shift and is 



2'-/ ' --a,-s\2'0 ^2 ^ 2"*^' 



(10) 



with 0's being "gauge" fields in terms of the fluctuations 
defined by {dt + avFd^)eaAx.t) = (jJaA^^i) ILs, 22]. To 
evaluate the effect of $o,s to the Josephson current, we 
have to average the Green's function in terms of auxiliary 
4> fields; g{x,t,t') =< g{x,t,t'\(f>) >0, where 
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(11) 



Here, 5ind = / dtdxLi is the action for the fluctuations 
induced by the interactions and forms a gaussian in terms 
of 9 fields. Considering that spin-charge separation is still 
valid, S'ind has the following form after unitary transfor- 
mation 



M = 











(12) 



where j = {/?, cr} specifies charge wave or spin wave fluc- 
tuations. One see that M corresponds to the inverse of 
Green's function of the fluctuations. Josephson current 
with s-wave symmetry is related only to charge bosons 
of the charge waves (i.e. Op^c = Yla s ^a,s/2) and phase 
bosons of the spin waves {9cr,p = ~Yl,a s'^^^asl'^) I3- 
Thus, it is sufficient to have the knowledge of the prop- 
agators of these components. The concrete forms in 
Fourier representation are 



MKp' - 2) 



+ Dl[k,uo) , (13) 



where 

D°{k,uj) = (vj-vf) 



""i^^S +D"Ak,^)], (14) 



Up' -I- VFVjk"^ 



(w^ — Vpk'^){uj'^ — v'^k'^) 



,(15) 
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FIG. 2: Josephson current in unit of Jc = evp/ L is potted as 
a function of phase difference for different KpS (Here, we set 
7 = 0.457r, L/Lt=QA, Ka=l). We assume gi\\ = gi±, thus 
R-^ = Uj. There are positive slopes at x = ^ indicating the 
existence of the vr-state. One can find critical currents are 
suppressed for repulsive interactions. 
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FIG. 3: Temperature dependence of in unit of Jc is plotted 
at 7 = O.SStt and — 1. For repulsive interactions, one can 
see that vr-state is strongly suppressed. Further, at Kp = 0.4, 
the generation of vr-state is entirely suppressed even at zero 
temperature. 



Here, Kj and Vj — UjVp are usual Luttinger parameter 
and propagation velocity of j component. 

Total Josephson current through the system can be ob- 
tained, summing up the contributions from all the ABSs. 
Because of the spin-charge separation, renormalization of 
the spin-up current and the spin-down one is exactly the 
same. Here we assume that excitation momentum unit 
of the charge and the spin bosons is infinitesimally small 
in LL connected with Fermi liquids. As a consequence, 
Josephson current for S-LL-S system has a following form 

J = —? 1— sm(TOx)cos(m7 — . , (16) 

L TT Lt smh(^) 

with Lt = vp I^'kT being thermal length. The influence 
of electron-electron interactions is included only in T, 
which has a following form 

T = (^^)Kr^+j^.-2 "p"^" ^^^'sinh^(2r7) 

^ Binh^"(2^I^)^i"^''"(ra)' 

(17) 

Here, D is a high energy cut-off; since we take into ac- 
count the energy smaller than superconducting energy 
gap, we assume D A ^^l- Eq.(16) shows that each 
AR process is renormalized by the fluctuations with en- 
ergies lower than vpjL This can be interpreted 
that the "emvironment" at the interfaces directly affect 
the AR. Note that the exponent of renormalization dif- 
fers from the one due to normal backscattering, e.g., at 
high barrier potential oc (7rT/D)^p"'+^- '"^ There- 
fore, we consider that it is Cooper pair tunneling pro- 
cess that the interactions renormalize. To see this explic- 
itly, we utilize conventional tunnel Hamiltonian Ht — 

Ea.a(^l*t(-f )^a,s(-f ) + Vl/t , ( f (| ) + H.C.) 



with \E'a_s(a;) and '0a,s(a;) = gi(afcF2;+e„,s) being annihi- 
lation operators in S and LL. At zero temperature, the 
second order perturbation of Ht at each interface, which 
corresponds to single AR process, is proportional to 

y < e'(''+T(-T'0)+«-i(-T.o)-e_i(i,t)-e+-r(#,t)) 

oc (^)^p"+^"-2 = T(r ^ 0). (18) 

This is in agreement with the result of Fazio et al. 0- 
One can verify that higher order contributions are con- 
sistent with Eq.(16) as far as the power law dependence 
is concerned jig. 

Current-phase relations for different values of Kp are 
shown in Fig 2 (Hereafter, we employ the parameters, 
such as A and L, modeling the experiment by Kasumov 
et al. 0). Owing to the Zeeman phase shift 7, the posi- 
tive slopes are found to appear at x = tt (vr-state). This 
is the same result as in S-F junction except the suppres- 
sion (enhancement) by the repulsive (attractive) inter- 
actions jl3|. However, the suppression by the repulsive 
interactions {Kp < 1) can play a key role as for the gen- 
eration of TT-state. We calculate the temperature depen- 
dence of the sub-critical current J,r which provides an 
indication of the existence of the 7r-state (see Fig. 3). In 
presence of the repulsive interactions, there is the regime 
for flnite 7 where the 7r-state does not occur even at ab- 
solute zero, although flnite Josephson current can flow. 
The condition for the existence of a 7r-state is given by 

00 

A - - ^ cos(m7)(-^)(^p"+^"-2)™' > . (19) 

m— 1 

In non-interacting case {Kp — K„ = 1), A = 0.5 for 
finite 7. In contrast, if the Coulomb interactions are 
strong enough {Kp < 0.2) as in CNTs 0, the metastable 
TT-state cannot be generated as such. 
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FIG. 4: The tunneling rate from the bottom of the vr-state 
is plotted for 7 — 0.35n (dotted line) and 7 — 0.457r (solid 
line). Declining line is the one plotted as a function of Kp at 
Kcr ~ 1 (left axis), while increasing one a function of Ka at 
Kp — 0.8 (right axis). 



This behavior can be verified also from Josephson po- 
tential Ej^XtJ)- As a criterion for determining the exis- 
tence of the metastable state, we estimate the transmis- 
sivity from the bottom of the 7r-state. T^r is numeri- 
cally calculated as a function of Kp and K^r within WKB 
method, and the result is shown in Fig. 4. One can see 
that there is a critical value of interaction strength, where 
reaches unity and the generation of 7r-state is prohib- 
ited completely. In the vicinity of the critical point, the 
tunneling probability becomes 

where C denotes junction capacitance. The disappear- 
ance of the metastable state is caused by the fact that 
the Zeeman splitting in ABSs becomes blurred due to the 
collective excitations arising from the interactions in LL. 
In practice, for repulsive interactions, the LDOS near the 
interfaces can be roughly estimated by 



n.s 



s\/D 



( 
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(e-e 
e — e 



(^)r(i-Cp-C.)cos(-{Cp + C.} 



(^^)^^)--r(i-.)cos(||c,.C. 



2^ 



(21) 



with the exponents Q = {Kj ^ + Kj — 2)/4: and rj = 
{K~^ + Kcr — 2)/4. Here S is the finite spectral weight 
of each ABS and r(a::) is Gamma function. We fur- 
ther fix the phase difference y between the superconduc- 
tors to zero for simplicity ,12j . The first term in right 
hand side show that the spectral weight is spread out 
through distribution function. In this light, the repul- 
sive interactions can be interpreted to raise the effective 
temperature of the junction system. In the second one, 
which originates in ^a,s of Eq.(lO), the repulsive inter- 
actions suppress the Andreev resonance. For attractive 
case {Kp > 1), on the other hand, it becomes easier 
to generate the metastable state. An attractive back- 
scattering, in general, cannot be neglected. However, it 
does not change the situation qualitatively, because the 
attractive back-scattering tends to enhance the super- 
conducting correlation [23 • Besides, in contrast to Kp, 
small K^ enhances the critical current, and makes the tt- 
state incident. This reveals a reflection of the symmetry 
of Cooper pairs in superconductors. 

Finally, we mention the effect of breakdown of spin- 
charge separation due to the Zeeman effect. We have fo- 



cused on the case where spin-charge separation is valid. 
This assumption may be acceptable for the system on 
electron doped GaAs which has small Lande's g-factor 

( 0.4) |||. If, e.g., InSb quantum wires (g^ 50 |2^ 'l 

are prepared, however, our approach needs some cor- 
rections. In such a case, the charge- and spin-density 
fluctuations exhibit coupled oscillation, whose coupling 
strength is proportional to r = {vpi^ — vfi)/{vf^ + ^Fi)- 
Therefore, M in Eq.(12) has off-diagonal elements. This 
leads to different renormalization for spin-up and spin- 
down currents, and the ratio of renormalization factor 
has the following expression 

where u+ denotes velocity renormalization of the normal 
mode the explicit form is given by 



2r') ± 



ulY + ir^iupKp + UcrKa^){upKp^ + UaK^) 



(23) 
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Since Eq.(22) does not include high energy cut-off D, 
/T| does not become so large. Thus the influence 
of the breakdown of the spin-charge separation is small, 
which may cause the circumstance that supercurrent is 
conveyed by singlet Cooper pairs. 

In conclusion, current-phase relations of S-LL-S 
junctions are analytically reexamined using functional 
bosonization. We have shown that the low energy excita- 
tions in LL can play crucial roles on a critical current and 
a generation of 7r-state if LL is adiabatically connected 
with superconducting reservoirs. We have also exam- 
ined the tunneling from the metastable 7r-state. When 
we regard the junctions as macroscopic two level sys- 
tems, the quantum leakage directly affects the stability 
of the metastable state. The decay rat e is given by F = 
2Be-^«-^'^('^'^)/2, where = dxy/CEj{x,Tr/2)/ie^ 
is the action of symmetric potential {B is the fluctua- 



tion determinant without zero mode 113, H^)- This in- 
dicates that the macroscopic nature is inevitably related 
with many body correlations in ID configuration. Our 
method can be easily applied to hybrid systems of super- 
conductors with other symmetries or interacting conduc- 
tor with higher dimensions I21L l29l | . Experimentally in 
semiconducting heterostructure, one can vary the carrier 
density and the effective interactions to some extent by 
gate voltage js^]. Further tunn el j unctions between ID- 
2DEG have been also examined jSJ . It may be feasible in 
the future to construct the system under consideration, 
using S-2DEG-S junction with a gate [3^ . 
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